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Regular quasiprobabilities are introduced for the aim of characterizing quantum correlations of 
multi-mode radiation fields. Negativities of these quantum-correlation quasiprobabilities are nec- 
essary and sufficient for any quantum correlation encoded in the multi-mode Glauber-Sudarshan 
P function. The strength of the method is demonstrated for a two-mode phase randomized squeezed- 
vacuum state. It has no entanglement, no quantum discord, a positive Wigner function, and a clas- 
sical reduced single-mode representation. Our method clearly visualizes the quantum correlations 
of this state. 
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Introduction. - Nowadays a variety of different no- 
tions of quantumness are discussed. It is of fundamental 
interest to classify the nature into one part showing clas- 
sical signatures, and another one dominated by quantum 
phenomena. In quantum optics, the most common no- 
tion of nonclassicality is established through the Glauber- 
Sudarshan P function [TJ [5] . If this distribution does not 
resemble a classical probability density, the correspond- 
ing state is called a nonclassical one [5]. 

Quantum entanglement is another notion of quantum- 
ness [U [5] . Although each individual subsystem can be 
quantum, for a separable state the correlations between 
the subsystems can be described by classical statistics [B] . 
Entanglement or inseparability can be fully characterized 
by negativities of optimized quasiprobability (QP) distri- 
butions [3 [8] , denoted as PEnt ■ 

In quantum information theory, a quantum correla- 
tion (QC) is also characterized by the so-called quantum 
discord [TU] . A non-zero discord describes the back- 
action of a measurement in one subsystem to another 
one. Discord is often considered as a general measure 
of QC beyond entanglement. Recently, however, it has 
been shown that the notions of QC based on quantum 
discord and on negativities of the P function are maxi- 
mally inequivalent 

The general characterization of the QCs of radiation 
fields has been formulated in terms of the space-time de- 
pendent P functional [H] . This yields a full hierarchy of 
inequalities for observable correlation functions. When- 
ever the described systems violate classical probability 
theory, the P functional is a QP representation of quan- 
tum light. In the following we restrict attention to equal- 
time measurements, so that the P functional simplifies to 
the ro-partite P function of the global quantum state p, 

p = J d 2n aP(a) \a)(a\, (1) 

where \a) — \a±, . . . , a n ) denotes multi-mode products of 
coherent states. One severe disadvantage of the P func- 
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tion consists in its strong singularities occurring for many 
quantum states, even in the single-mode case. As a con- 
sequence, in general this function is experimentally not 
accessible and hence of limited practical value. Only in 
special cases the P function is regular [13] and can be 
measured |14) . 

The concept of phase space representations has been 
generalized to the s-parametrized QPs [TS]. The latter 
include popular examples: the P function (s = 1), the 
Wigner function (s = 0), and the Husimi Q function 
(s = — 1). Diminishing the s parameter, the QPs become 
more regular, but less sensitive for testing quantumness. 
For squeezed states, all s-parametrized QPs are either 
positive or irregular. Especially the Wigner function is 
popular, since it is easily obtained in experiments, e.g. 
for quantum light, molecules, and trapped atoms [TBl - 
20 . Further generalizations of the QP methods were in- 
troduced in [21] . 

For the single-mode case, nonclassicality QPs, PncI, 
have been introduced [22] . They are regularized versions 
of the highly singular P functions. For any nonclassi- 
cal single-mode state, they show negativities and can be 
directly obtained from experimental data [22 GH [38] . 
Hence Pn c i is a powerful tool for the full experimental 
characterization quantum effects of single-mode fields. 

In the present contribution, we generalize the regular- 
ization of the P function for multi-mode light. The re- 
sulting QP distribution, Pqc, uncovers any QC occurring 
in the multi-mode P function. The former can be directly 
sampled from experimental data obtained by multi-mode 
homodyne detection. The method is applied to a phase- 
randomized two-mode squeezed vacuum, which is not en- 
tangled, has zero quantum discord, a positive Wigner 
function, and classical reduced subsystems. Even in this 
case, our method uncovers QCs contained in the highly 
singular P function, via negativities of the bipartite fil- 
tered QP distribution Pqc- The QCs are clearly visible, 
even when the other signatures of QCs do not persist. 

Regularized multipartite Pqc function. - To be useful 
for experiments, quantumness criteria must be based on 
well-behaved functions. Hence a regularization of the 
highly singular Glauber-Sudarshan P representation is 
required. We can construct a multi-mode regular Pqc 
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function in the form of a convolution 

P QC (a;w)= J d 2n a' P(a-a')a w {a'), (2) 

where Q w is a suitable function, which we are going to 
construct. The occurring width parameter w provides 
the property Pqc(cx]w) — > P(<x) for w —> oo. A direct 
sampling formula for measured quadrature data will be 
given. 

We start from the multi-mode characteristic function. 
The n-mode characteristic function $ is defined as the 
Fourier transform of the P function 



$(/3) 



(3) 



The characteristic function $ is always a continuous func- 
tion, independently of singularities in P(a). Single-mode 
nonclassicality criteria based on the characteristic func- 
tion have been introduced [57] and applied [2"5H50] . 
Note that the n-mode characteristic function is bounded 
in the form |$(/3)| < exp(|/3| 2 /2). 

The convolution in Eq. ^ is a point-wise product in 
Fourier space. Hence, we consider the so-called multi- 
mode filtered characteristic function $q C (/3; w), 



(4) 



k=l 



being the Fourier transformed of Pqc for < w < oo. 
Our choice for the filter Cl(/3) is the autocorrelation func- 
tion 
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Finally, we get the regularization kernel in Eq. ([2| by the 
inverse Fourier transform as 



m«) = n (~ /d ; 

k=l v 



p fl(P/w) e 
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For the single-mode case, such a filter Q(P/w) in 
Eq. ([5]) has been introduced and characterized in 
Ref. [22] . It has been shown that this filter suppresses the 
exponentially rising behavior of exp(|/3| 2 /2). In addition, 
it has been also shown that this kind of filter belongs to 
the class of invertible filters, Cl(P/w) 7^ 0. Hence, decon- 
volution of Eq. |2]) yields the P function, and thus the 
full state p, from the QP distribution Pqc. 

Let us comment here on the structure of the regular- 
izing function Cl w defined in Eq. This multi-mode 
function can be written as a product of single-mode func- 
tions which do not depend on the state, the proof is given 
in the appendix A. The regularizing function must not in- 
troduce any quantum correlation to Pqc , which is absent 
in the Glauber-Sudarshan P function. For any nonclas- 
sical state which includes some QC, there exist values w 



and a for which Pqc (a; w) < 0. The other way around, 
a state is classical, if for all values uu and a the func- 
tion Pqc (a; w) represents a classical probability density. 
This means that we may identify QCs via uncorrclated 
filtering. 

Sampling of Pqc- - For applications of our method, 
a direct sampling formula is useful which yields Pqc from 
the measured data. We consider multi-mode homodyne 
detection, which gives a sample of N measured quadra- 
ture values {(x[j} 7 (p[j})}f =1 , with (x k [j],<f k [j]) being the 
j-th measurement outcome in the fc-th mode. Since our 
multi-mode regularizing function Q w can be written as a 
product, we obtain the sampling formula 

N n 

P QC {a;w) = — II f( x k[j]^k[j],a k ;w). (7) 
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The needed single-mode pattern function reads as 
-1 p+00 

f(x,ip,a;w) =— dbbe b2/2 n(b/w)x 



(8) 
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it has to be calculated only once. 

In order to provide the confidence of the negativity of 
the sampling in Eq. Q , we need to estimate the statisti- 
cal uncertainties. The empirical variance of the sampling, 

[Pqc(<*; w)}, is derived in the appendix B. The statis- 
tical uncertainty of the sampled estimate in Eq. Q is 



A[P Q c(a; 



a [Pqc (a; w)} 



N 



(9) 



Thus one can obtain the confidence of a negativity as 



C(ac; w) 



\P QC (a;w)\ 
A[P Q c{a;w)} 



(10) 



for a given a with Pqc (a; it;) < 0. One may choose the 
value w > 0, such that C(a\ w) becomes maximal. 

Beyond entanglement and discord. To show the 
strength of our method, we consider a realistic exper- 
imental scenario for generating a bipartite continuous 
variable state. The inputs of a 50:50 beam splitter are 
equally squeezed states in orthogonal quadratures, c.f. 
Fig. [T] In the output ports of this setup one obtains a 
two- mode squeezed vacuum state (TMSV). Such states 
are entangled and have non-zero quantum discord. 

The situation is drastically changed, if phase random- 
ization, indicated by dep in Fig. [T] occurs in one of the 
output ports. The entanglement properties of the state 
depend sensitively on the dephasing [5] . In the case of an 
equally distributed phase, Sip — 2ir, the resulting output 
state is a fully phase randomized TMSV state, 
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FIG. 1: Experimental setup for the generation of a phase ran- 
domized two-mode squeezed vacuum. A 50:50 beam splitter 
combines squeezed-vacuum states, which are squeezed in or- 
thogonal quadratures. One of the output channels is fully 
phase randomized, Sip = 2n. 



with < p < 1, p being related to the squeezing param- 
eter of the initial input fields. 

For the following it is important to consider properties 
of this state which are closely related to its QCs. First, 
the phase randomized TM SV state is a convex mixture of 
the product states <g> \n)(n\. Hence, it is classical 

with respect to the property of entanglement. Second, 
due to the orthogonality of the Fock states, (n\n') = for 
n 7^ n', this state has no QCs in the sense of the quantum 
discord. So far our state belongs to the class of quantum 
states considered in [TT] with the aim to demonstrate the 
inequivalence of QCs based on the P function and the 
quantum discord. In our context, we require the following 
additional properties. 

Third, let us consider the reduced density operator, 
/5 ro d. Due to the symmetry of this state with respect to 
the interchange of both subsystems, we find that 



p rod = Tr ap = Tr B p = - p)p n \n}(n\ 



71=0 



(12) 



This state is a thermal one with a mean photon number 
fhh = p/(l — p)- Therefore it shows a classical behav- 
ior with respect to the reduced single-mode states. This 
property insures that any identified signature of quan- 
tumness exposes a true QC effect. Fourth, the full state 
p is a mixture - due to the phase randomization - of 
two-mode Gaussian states. Hence it has a non-negative 
Wigner function, which is shown in Fig. [2j This QP 
distribution, which was often applied in experiments |16l - 
[20] , is improper to visualize QCs of this state by attaining 
negative values. 

Summarizing the features of the fully phase random- 
ized TMSV, we have introduced an experimentally real- 
izable quantum state with the following properties: 

1. no entanglement; 

2. zero quantum discord; 




FIG. 2: The Wigner function is shown for a fully phase ran- 
domized two-mode squeezed vacuum state, with p — 0.8. It 
is non-negative in the full two-mode phase space. 



3. classical reduced single-mode states; 

4. non-negative Wigner function. 

Despite these strong signatures of classicality, this state 
will be proven to show clear QC effects. 
The P function of our state is given by 



P(a A ,a B ) = J2(l-P)P n Y, (")(?)< 11 



n=0 



k,l=0 



k+l 



x m d*J k a *J(a A )dl B d l a *J(<XB). (13) 

It is a non-Gaussian state that can be experimentally 
generated as outlined above. To prove that this state de- 
scribes QCs, we have to visualize the negativities of this 
strongly singular P function. If this attempt succeeds, 
we have shown that there are QC between the subsys- 
tems A and B going beyond entanglement and quantum 
discord. 

For the proof of the existence of QCs, the regularized 
two-mode QP, Pqc, is calculated from Eq. (2]) and (13) 
together with Eq. ^ . The result in Fig. [3] clearly shows 
that Pqc becomes negative for properly chosen argu- 
ments. In view of the properties of the state listed above, 
these negativities are a direct proof of the quantum na- 
ture of the correlations between the two subsystems A 
and B. Although the Wigner function in Fig. [2] contains 
the full information on the quantum state, it does not 
directly visualize the QC properties of the state under 
study. The regular QP distribution Pqc, on the other 
hand, directly displays the QCs within the considered 
randomized TMSV state by attaining negative values. 

Here the question may arise, whether the quantum cor- 
relations of the fully phase randomized TMSV state can 
be used in any application. It opens possibilities for the 
heralded control of the number of photons in an arbitrar- 
ily chosen time interval. More generally, even the time 
sequence of a train of photons could be controlled by the 
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FIG. 3: The regularized quasiprobability Pqc is shown for 
a fully phase randomized two-mode squeezed vacuum state, 
for p = 0.8 and w = 1.5. Clear negativities visualize the 
nonclassical correlations of the state. 

quantum correlations inherent in such a state. Presently, 
the control of photon states for larger number of photons 
is a demanding problem. 

Conclusions.- A regularized multimode version of 
the Glauber-Sudarshan P function is introduced: the 
quantum-correlation quasiprobability Pqc, which is pa- 



rameterized only by a single width parameter w. 
Through its negativities, it directly visualizes quantum 
correlations included in any multimode quantum state. 
The negativities occur, if and only if the state has a 
nonclassical multimode Glauber-Sudarshan P function. 
The method of regularization is universal, since it does 
not depend on the considered state. Via direct sam- 
pling of experimental data by multimode balanced homo- 
dyning, one can directly access the quantum-correlation 
quasiprobabilities . 

The method is applied to a fully phase-randomized 
two-mode squeezed-vacuum state. This state is shown 
to be classical with respect to other established no- 
tions of quantum correlations and it is locally classical 
in each mode. However, the negativities in the bipar- 
tite quantum-correlation quasiprobability clearly uncover 
the quantum correlation properties of this state, beyond 
the existence of quantum entanglement and quantum dis- 
cord. 
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Appendix A: Regularization 

Here, we prove the properties of the multi-mode filtering procedure. For the single-mode scenario, we refer to the 
appendix in Ref. [32]. First, we show that the a product filter is sufficient to identify any kind of nonclassicality 
contained in the multi-mode P function. Second, we show that the filtered quasiprobability Pqc is a smooth function. 

1. Product filters 

For convenience, we will argue in the Fourier domain. The multi-mode filtered characteristic function (CF) reads 

as 

n 

$ QC (/3; u;) = $(/3) jj £l(f} m /w), (Al) 

m— 1 

with Q being the continuous function 

(1(0) = {^j' 2 J d 2 /?'e-l^'lVl^\ (A2) 

with lim^-^oo Vl((3/w) = 1. 

Bochner's theorem can be given in the following form [26l [35]. The CF 3>qc is the Fourier transformed of a 
probability distribution, iff for all finite sequences £ C and all finite sequences (3 k € C" holds that 

E K f*> $ QC 09* - 5 «>) > 0: (A3) 

together with the normalization <I>qc(0;u>) = 1 and the symmetry $qc(/3;u')* = $qc( — We may rewrite the 
matrix (<I>Qc(/3 fc — f3 k ,; w))u,h' in the following form: 

($Qc(/3 fe - P k r,w))k,k> = (®(Pk-Pk>))k,k'°to(\Pi t k-Pi,k']/ w h t k> ° ■ ■■ ° tt{[Pn,k - Pn,k>]/ W )k,k>, (A4) 
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where o denotes the entrywise or Schur product. It is well-known that the Schur product of positive definite matrices 
is again positive semidefinite, and the matrix (l)k,k' is the identity element of the Schur product. 

Let us assume that our state is a classically correlated one. This means that the multi-mode P function is a non- 
negative probability density, which implies that $ fulfills Bochner's theorem, cf. Eq. (A3), or simply {${(3 k — /3k>))k,k' 
is a positive semidefinite matrix. From the definition of f2, see Eq. (A2), it is clear that Q{[ftm,,k~ Pm,k']/w)k,k' (f° r the 
vector components m = 1, . . . , n) is also positive semidefinite. In the case w — > oo, we have Q([/3 TO> fc ~ /?m,fc']/w)fc,fc' — > 
(l)fc,fc'- This means that for any w > 0, we have a positive semidefinite matrix ($ qc(/3 fc — (3 k ,; w))k k'- 

Let us now consider a quantum-correlated state. Thus, the condition in Eq. (A3) cannot be true. This means 
that there exists a finite sequences € C and a finite sequences j3 k G C", such that the quadratic form ($Qc(/3 fe — 
f3 k ,;w))k,k' is not positive semidefinite. The continuity of fl for w — > oo implies f2([/3 m ,fc — Pm,k']/w)k,k' ~^ (l)fc,fc' ; we 
get that there exists a finite w > such that {^QciPk ~~ @k''i w ))k,k' is not a positive semidefinite matrix. Thus, the 
inverse Fourier transformed Pqc(c*; w) cannot be classical. This holds true for an interval of width parameters from 
the particular w to oo. 

Concluding this first part, we have seen that the product filter for a particular choice of w exhibits any kind of 
quantum correlations within the state. The other way around, a state is classical, iff for any w > 0, we have a 
non- negative -Pqc- 



2. Smoothness 



It is needed to consider the complex plane, given by its real and imaginary parts, as a two-dimensional real space, 
C = Ixl. Our aim is the proof that Pqc( a ', w ) is a smooth function, Pqc € C°°(C™, K). In Ref. [22, it was shown 
for the single-mode case (n — 1) that the filtered quasiprobability is a square-integrable function. Here, we want to 
go further and prove that the multi-mode filtered function is smooth of any w > 0. 

The properties of the characteristic function are basically given by its definition as $(/3) = (:D(f3):) , where :D((3): = 
e p a e -/3 a re p resen t s th e multi-mode normally ordered displacement operator, see e.g. [3B]. The relation to the 
unitary displacement operator is given by e~'^' / 2 :l)(/3): = D{j3). In addition, we have a continuity in the form 
D{/3 + 0) = c-' ilnl ^'^^D{0)D(/3'). Thus, we have that $(/3) is continuous and bounded in the form 

|$(/3)| <e l/3|2/2 . (A5) 

In the following, we need to prove that the Fourier transformed of 'I'qc is smooth. In addition, we may summarize 
some properties of the Fourier transform J 7 . The Fourier transform is an automorphism of square-integrable functions 
L 2 (C",C). The inverse Fourier transform of an auto-correlation function is J-Aut[g] = |.F<7| 2 . It is also important to 
recall the Sobolev's lemma, cf. e.g. [37]: If |/3| fc <?(/3) € L 1 (C TI ,C) for k £ N, it follows that the Fourier transform is 
continuously differentiable, Tg £ C fc (C n ,C). 

Theorem 1. Let w(/3) be a positive function, which satisfies w(/3)e"" 3 ' € L 2 (C n ,R) for all u > 0. Then Pn(a; w) e 
C°°(C n ,R) for all w > 0. 

Proof. We have 

P QC (a; w) = J-[$(/3)fi(/3/ w )](a). (A6) 

It is sufficient to prove that $ n (/3;w) = $(/3)fi(/3/w) G L 2 (C",C), and |^| 2fc $(/3)fi(/3/w) e V-(C n ,C) for all w > 
and k € N. We have 

||/3| 2fc $ QC (/3;^)| <C u \(3\ 2k e^ 2 / 2 e^ 2 / w2 (A7) 
<C„fc!e( 1+ 5-^)l/3| 2 . (A8) 

We can choose a u > 3w 2 /2, such that the left-hand side is bounded by a rapidly decaying function, which is in 

i 1 (c™,c)n J L 2 (c™,c). ' □ 

Note that the considered u)(j3) — e - '' 3 ' fulfills the requirements. We can bound it as o;(/3)e u '^' = e u / 4 e C I 1 -«/ 2 ) 



Appendix B: Sampling error estimation 

In the following we present the calculation of the empirical variance of the sampling. The error estimation for 
sampling methods of one-mode systems could be founded in |38j and the references therein. The statistical uncer- 
tainties below calculated allow us to determine the confidence of any value of the filtered multi-mode quasiprobability 
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distribution Pqo The variance of the PQc(ot;w) function is 

a 2 [P QC (a;w)} = E [P QC (a;w) 2 ] -E[P QC (a;w)] 2 . 
The sampling formula for Pqc(ck; w) is given in terms of a product of the so called pattern functions as 



N n 



Pqc(a;w) = — ^ J| f(xk\j], Vk\j], a k ; w), 



(Bl) 



(B2) 



j=i k=i 



where the latter are constructed based on a sample of N measured quadrature values {(x[j], <p[j])}j =1 , with 
(xk[j],fk[j}) being the j-th measurement outcome in the fc-th mode. 
The expectation value of the second moment of Pqc(ol;w) is 



N 



j=l k=l 



1 



E[P QC (a;w) 2 ] =E 
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\ J f{xk[i]^k[j],ak]w) 



I iV n 

V X II f(, x k\j],Vk\j],akV> 



N n 



j'=ifc=i 



N 



.k=i 



AT 



.fc=i 



E 



II /(^feb"']) Vfe b'W; 
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,fe=l 
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iV2 
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|1 f 2 (xk,(Pk,a k ;w) 



k=l 



+ N(N- 1)E 



II f(xk,<Pk,a k ;w) 



Lk=l 



(B3) 



We just decompose the two original sums into all contributions with j = f plus j ^ j' and took into account the fact 
that individual measured points {{x[j},(p[j})} are stochastically independent and equally distributed. The expectation 
value of the regularized P function is 



1 N 

E[P QC (a ;U ;)]=-^E 



3 = 1 



f(xk[j],fk[j},a k ;w) 



k=l 



1 

N 



NE 



JJ f(xk,<Pk,a k ;w) 
Lfe=i 



(B4) 



Finally the empirical variance is given in terms of the expectation value of the product of pattern functions and the 
product of its squared value as 



a 2 [P QC (a;w)} 



1 

N 



E 



{[ f 2 (x k ,<f k ,a k ;w) 



Lk=l 



E 



I 2 



n f( x k^k,a k ;w) 



Lk=l 



(B5) 



For the particular case of a two-mode system this variance looks like 

<T 2 [PQ C (a;w)] = |e [/ 2 (xi, (^i,ai; w) / 2 (x 2 , (p 2 ,a 2 ;w)] - E {f(x 1 ,Lp 1 ,a 1 ;w) f(x 2 , <fi 2 , u 2 ; u>)] 2 } , (B6) 

and for our particular example of the fully phase randomized two-mode squeezed vacuum state the phase-independent 
pattern function is 

f(x,a;w) = - dbbe b ^ 2 e tbx Jo{2b\ ai \)n{b/w). (B7) 
^ Jo 

Then the second moment and the expectation value of the regularized P function for this state are respectively 

(n \ 4 /.00 />oo 
-J ^ d&i6iJ (26i|ai|)fi(&i/t«) y dftifti J (26 / 1 |a 1 |)O(6' 1 / U ;) (B8) 

roc poo 

/ d& 2 6 2 J (2b 2 \a 2 \)n(b 2 /w) / d6 2 6' 2 J (2fe 2 |a 2 |)r!(6 2 / W )e bl6 'i+ b2b 2$(6 1 -6' 1 ,6 2 -6y 
Jo Jo 



and 



(rj \ 2 /.00 />00 
-J y d66J (26|ai|)n(6/«;) y db'b' J (2b'\a 2 \)fl(b' /w)^{b,b'). (B9) 
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The last expression is just Pqc [cx-i , Q-i \ w). The variance written in terms of ( B8 ) and ( B9 ) is applied in the calculation 
of the confidence. I 
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